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Massive, luminous stars reaching the Eddington limit in their interiors develop very dilute, 
extended envelopes. This effect is called envelope inflation. If the progenitors of Type Ib/c 
supernovae, which are believed to be Wolf-Rayet (WR) stars, have inflated envelopes then 
the shock breakout signals diffuse in them and can extend their rise times significantly. We 
show that our inflated, hydrogen-free, WR stellar models with a radius of ~ R© can have shock 
breakout signals longer than ~ 60 s. The puzzlingly long shock breakout signal observed in the 
Type Ib SN 2008D can be explained by an inflated progenitor envelope, and more such events 
might argue in favour of existence of inflated envelopes in general. 


1 Stellar envelope inflation 

Inflation refers to the extremely dilute, loosely- 
bound envelopes that massive, luminous stars de¬ 
velop in the course of their evolution (Ishii et al. 
1999; Petrovic et al. 2006; Sanyal et al. 2015, Langer 
et ah, this volume). A typical inflated envelope is 
shown in Fig. 1, which corresponds to an evolved 
hydrogen-free 7.8 M© model burning carbon in the 
core. Such ubiquitous hydrostatic structures result 
from stellar layers reaching the Eddington limit in 
the interior of the star, facilitated by the peak in 
opacity at T « 200 000 K, the so-called iron opacity 
bump. For more details on the appropriate defini¬ 
tion of the Eddington limit, we refer to Sanyal et al. 
(2015) and Langer et al. (this volume). 



Fig. 1: Density profile of a 7.8 M© hydrogen-free stel¬ 
lar model with Teff = 82 700 K and log(L/L 0 ) = 5.35, 
showing an inflated envelope (shaded blue) and a density 
inversion. The radial co-ordinate has been normalised by 
the non-inflated core radius (rcore) of 0.6 R© . 

Stellar layers approaching the Eddington lumi¬ 
nosity become convectively unstable (Langer 1997) 
but convection, of course, does not contribute to 
the radiative acceleration. However, the fraction 
of the total flux transported by convection in these 


WR stars, according to standard MLT, is negligibly 
small. Hence, whether or not the layers will hit the 
Eddington limit and inflate the envelope, does not 
depends critically on the theory of convective energy 
transport used to model the stars. 

As we can see from Fig. 1, the mass contained in 
the inflated envelope is a tiny fraction of the stellar 
mass, AMenv ~ 10“® M© . However, for cool super¬ 
giants, the mass in the inflated envelopes can be as 
high as a few solar masses, see Sanyal et al. (2015) 
for details. 


2 Stellar models of WR stars 

We modelled WR stars as non-rotating, hydrogen- 
free helium stars and evolved them through core¬ 
helium buring to the onset of core-carbon burning, 
using the 1-D stellar evolution code BEC (see Heger 
et al. 2000; Brott et al. 2011; Moriya et al. 2015, 
for details). We computed two sets of models, with 
initial masses 10 M© and 12 M© , and a metal com¬ 
position that of the Sun (Brott et al. 2011). In par¬ 
ticular, we used the OPAL opacity tables (Iglesias 
& Rogers 1996) for our computations and used the 
standard Mixing-Length-Theory (MLT) to model 
convective zones, with omlt = 1-5. The mass-loss 
rate prescription from Nugis & Lamers (2000) was 
applied. The resulting structure (e.g. Fig. 1) will 
resemble a SN Type Ib progenitor because no signif¬ 
icant changes in the stellar structure is expected in 
the remaining lifetime until the explosion. 

The evolutionary tracks of the aforementioned 
models, in the H-R diagram, is shown in Fig. 2. 
Both the helium zero age main sequence models 
are marginally inflated. As a consequence of the 
applied mass loss, they decrease in luminosity dur¬ 
ing the core helium-burning phase. After helium is 
exhausted in the core, the stars begin to contract 
and become hotter and brighter. The models are 
hardly inflated during this phase because the Fe- 
group opacity bump is only partially contained in¬ 
side the stars. These are probably related to the WO 
stars (cf. Tramper et ah, this volume). When the 
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helium-shell ignites, the evolutionary tracks even¬ 
tually move toward cooler temperatures because of 
the mirror principle. As the models become cooler, 
they become significantly inflated and exhibit a pro¬ 
nounced core-halo structure, as shown in Fig. 3 (cf. 
Moriya et al. 2015). 

3 Shock breakout in inflated 
stellar envelopes 

The shock breakout occurs when the dynamical 
timescale of the shock propagation in the unshocked 
envelope, i.e. tdyn — Ai?/ush becomes compara¬ 
ble to the diffusion timescale in the envelope, i.e. 
tdiff — tAR/c, where Vsh is the shock velocity and r 
is the optical depth in the remaining unshocked enve¬ 
lope (Weaver 1976). Therefore, the shock breakout 
condition can be expressed as: t ~ c/vsh- 



Teff[kKl 


Fig. 2 : Stellar evolutionary tracks of our 10 Mq and 
12 Mq solar-metallicity helium-star models from the He- 
ZAMS until core-carbon ignition. 


For a model in which the shock breakout condi¬ 
tion is fulfilled very close to the stellar surface, AR 
is very small, such that tdiff ^ ^ic, where tic — R*/c 
is the light-crossing time. Assuming a shock velocity 
of 20 000 km s“^, in the case of the polytropic model 
shown in Fig. 3, the light-crossing time tic is 5.6 s 
whereas the diffusion time is only 0.5 s, meaning that 
the shock breakout duration is determined by tic. 

On the other hand, in an inflated stellar model, 
the shock breakout happens relatively deep inside 
the envelope, making AR quite large, such that 
tdiff tic. In Fig. 4, we plot the optical depth in 
the interior of our stellar models and show that AR 
can be as large as ~ Rq . In such a situation, the rise 
time of the shock breakout signals will be determined 
by the diffusion time. The subsequent light curve is 
expected to decline exponentially with an e-folding 
timescale of the diffusion time due to the photon dif¬ 


fusion in the shocked envelope (Moriya et al. 2015). 
Table 1 summarises the different timescales involved 
in inflated and un-inflated models. 



Fig. 3: Density structure of our last computed models, 
with initial masses of 10 Mq and 12 Mq . The final 
masses are indicated in parentheses, in units of a solar 
mass. The density structure of a 7 Mq polytropic model 
(polytropic index = 3) without an inflated envelope is 
also shown. 


SN 2008D: The X-ray light curve of the Type Ib 
SN 2008D was serendipitously observed by Soder- 
berg et al. (2008) with the SWIFT satellite, and 
is the only direct detection of a SN shock break¬ 
out so far. The total shock breakout duration was 
300 seconds with a rise time of ~ 60 seconds. If we 
assume the shock brekout duration was dominated 
by tic, as is done in the literature, then it implies 
a WR progenitor radius of 130 R© , which is much 
larger than predicted by stellar evolutionary mod¬ 
els. On the other hand, our inflated model with a 
radius of ^ 2 Rq can naturally explain the long rise 
time of SN 2008D. The estimated SN ejecta mass 
is 3 — 7 Mq (Soderberg et al. 2008; Mazzali et al. 
2008; Bersten et al. 2013) which is consistent with 
the final masses of our models (after subtracting the 
remnant mass of the neutron star). The mass-loss 
rates of our models (7 x 10“® Mq yr“^) are also con¬ 
sistent with that estimated from radio observations 
(7 X 10“®MQyr“^, Soderberg et al. 2008). 

The shock velocity at shock breakout has been 
suggested to be ^ 30000kms“^ for WR stars (e.g., 
Nakar & Sari 2010), but these studies assume a 
steeply declining density profile as in the case of the 
polytropic model. Hydrodynamic modelling of shock 
propagation in inflated envelopes is needed to esti¬ 
mate the shock velocity at breakout. 

There have been other models proposed for ex¬ 
plaining the long shock breakout signal of SN 2008D, 
like the thick wind model by Svirski & Nakar 
(2014) which assumes a high mass-loss rate of ^ 
10“‘^MQyr“^ a few days prior to the explosion. 
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and the supernova ejecta expanding through the 
optically-thick wind. 

Tab. 1: A comparison of the properties and the relevant 
timescales of our computed models. The shock breakout 
rise time of SN 2008D was ~ 60 seconds. 


Mass ( Mq ) 

tic (s) 

^diff (s) 

Rn (R©) 

6.8 

5.67 

14 

2.43 

7.8 

5.38 

50 

2.31 

7.0 (Polytropic) 

5.60 

0.5 

2.40 



Radius [R; ] 

Fig. 4: The run of optical depth r in our models. For 
computing r in the polytropic model, a constant opacity 
coefficient of 0.2cm^g“^ was used. The green dashed 
line at r = 15 is drawn for convenience. The blue and 
red dotted lines show where in the stellar envelope shock 
breakout happens. 


4 Conclusions 

We have shown that if WR type supernova progen¬ 
itors have inflated envelopes, then the rise times of 
supernova shock breakout signals can be extended 
because the shock breakout can then occur within 
the low-density envelopes. The long diffusion time 
of the inflated envelopes makes the shock breakout 
rise times long. Even if a SN progenitor has a ra¬ 
dius on the order of the solar radius whose light¬ 
crossing time is a few seconds, the rise time of the 
shock breakout signals can be ~ 60 s because of the 
inflated envelope. Our inflated model can simulta¬ 
neously explain the mysterious long shock breakout 


rise time, the ejecta mass, the mass-loss rate esti¬ 
mate right before explosion, and the shock velocity 
in SN 2008D. 

The inflated envelope is a generic feature of lumi¬ 
nous, mass stars and can have an array of observa¬ 
tional consequences. It leads to cooler effective tem¬ 
peratures, higher spin-down rates in rotating stars, 
and less massive supernova progenitors. It has been 
argued in the literature (Grafener et al. 2012) that 
inflated WR envelopes can address the WR radius 
problem. A recent study of galactic OB stars by Cas¬ 
tro et al. (2014) indicates that inflated envelopes do 
exist in nature. Furthermore, core-hydrogen burn¬ 
ing stars evolving to temperatures below ^ 8000 K 
can be expected to have massive inflated envelopes 
and might be related to LBV giant eruptions (Sanyal 
et al. 2015). 
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